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Abstract. Charged particles of extensive air sliow- 
ers (EAS), mainly electrons and positrons, initiate 
the emission of fluorescence light in the Earth's 
atmosphere. This light provides a calorimetric mea- 
surement of the energy of cosmic rays. For recon- 
structing the primary energy from an observed light 
track of an EAS, the fluorescence yield in air has to 
be known in dependence on atmospheric conditions, 
like air temperature, pressure, and humidity. Several 
experiments on fluorescence emission have published 
various sets of data covering different parts of the 
dependence of the fluorescence yield on atmospheric 
conditions. 

Using a compUation of pubUshed measurements, 
a calculation of the fluorescence yield in dependence 
on altitude is presented. The fluorescence calculation 
is appUed to simulated air showers and different 
atmospheric profiles to estimate the influence of the 
atmospheric conditions on the reconstructed shower 
parameters. 
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emission, temperature-dependent colUsional cross 
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I. Introduction 

The number of emitted fluorescence photons at the air 
shower can be written as 



r(A,P,T,e) 



(1) 



dXdX ^ ' ' ' ' dX ' 
where Y(X, P,T,e) is the fluorescence yield in depen- 
dence on wavelength A, air pressure P, air temperature 
T, and vapour pressure e. The deposited energy of the 
secondary particles is denoted as dE^^^^/dX . 

In the last couple of years, a lot of effort has been 
put on the investigation of atmospheric dependences on 
nitrogen fluorescence in air |[T]. The fluorescence yield 
Y\ can be written as 

where $° is the fluorescence efficiency at zero pressure, 
P is the air pressure, and P' is the characteristic pressure 
for which the probability of coUisional quenching equals 
that of radiative de-excitation. The index v indicates the 
vibrational level of the exited state. Several groups have 
already investigated aspects of the fluorescence emission 
from nitrogen molecules in air (e.g. Bunner 121, David- 
son & O'Neil la, Kakimoto et al. |4|, MACFLY l5i 



and FLASH |6|). In addition there are various ongoing 
experimental activities, e.g. AIRFLY E), JS), ID, (TOl, 
Nagano & Sakaki et al. ITT], Gil, AirLight IH and 
Ulrich & Morozov et al. (14^. One major goal of all 
experiments is to obtain an absolute fluorescence yield 
y^* = • A//ic either for the main contributing band 
at 337.1 nm or for the entire spectrum in the range 
of interest between about 300 - 420 nm. repre- 
sents the intrinsic radiative de-excitation of the nitrogen 
molecules. However, in gas like air quenching processes 
have to be taken into account because the rate of 
radiative de-excitations is reduced by collisions between 
excited nitrogen molecules and further molecules in the 
gas. These quenching processes depend on atmospheric 
conditions and are described by (l+P/P^)^^ in Eq. (|2]l. 
Accounting all currently known effects, we can write 
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with To.i; as the mean life time of the radiative transition 
to any lower state, the index v indicates again the 
vibrational level of the exited state as for P^, Na is 
Avogadro's number, R is the universal gas constant. 
Pair is the air temperature, k is the Boltzmann constant, 
Cvoi is the fractional part per volume of the relevant 
gas constituents, and is the mass per mole where 
X stands for the relevant gas constituents. Up to now, 
the collisional cross sections cn^.d have been taken 
as temperature-independent even though it was known 
from theory that there has to be a temperature depen- 
dence. Recently, first experiments could confirm this 
dependence for nitrogen-nitrogen and nitrogen-oxygen 
quenching. The temperature-dependence of the nitrogen- 
vapour quenching has not been measured yet. First 
estimates indicate only minor importance with an effect 
of less than 1% change in the reconstructed energy of 
an air shower lITSl . An independent measurement of the 
temperature-dependent collisional cross sections in air 
has been performed quite recently. First analyses of data 
indicate compatible results with the measurements from 
AIRFLY and will be published soon ifTsl . 
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Adopting this description of fluorescence emission for 
air shower reconstruction, we have to apply atmospheric 
profiles for temperature, pressure, and vapour pressure. 
This cannot be provided by simple atmospheric models 
as these usually do not include vapour profiles. However, 
profiles obtained with meteorological radio soundings do 
provide all necessary quantities liT6l . 

II. Fluorescence Models in Reconstruction 

For this study, we could use the simulation and 
reconstruction framework Off line ifTTll of the Pierre 
Auger Observatory fTF]. Within this framework, we 
could obtain standard monthly models for the area of 
that observatory which do not include water vapour 
profiles |16|. Additionally, we had access to 109 actual 
nightly atmospheric profiles from local radio soundings 
that cover all conditions within a year. One of the 
advantages of the framework is that it features many 
implementations of different fluorescence models which 
can easily be interchanged. 

The first implementation of a fluorescence model 
in Off line referred to as K96, is based on measure- 
ments by Kakimoto et al. |4|. The fluorescence yield 
is parametrised in dependence on deposited energy 
and on altitude by considering the pressure and VT- 
dependences. The second fluorescence model, N04, 
has the same functional form of parametrisation and 
describes data from Nagano et al. [191 , LHJ- These 
measurements provide spectrally resolved data for 15 
wavelengths between 300 and 430 nm. Also in this 
description, only the pressure and \/T-dependences are 
considered. The third fluorescence description in Of fline 
is given by the AIRFLY Collaboration in 2007, labelled 
with A07. The fluorescence yield is given as 19J 

1 + _S! 

P'(A,To)^T/To 

^Po^To fluorescence yield at 337.1 nm as measured 
at their standard experimental conditions which are Pq 
= 800 hPa and To = 293 K. The other transitions 
have been measured relatively to that at 337.1 nm and 
are given by /p^ j,^. Overall, 34 transitions could be 
resolved between 295 and 430 nm. Since the absolute 
calibration of this experiment is still under study, ¥337 is 
normalised to the corresponding value of N04. It should 
be pointed out that the description in this model can 
easily be expanded to account for vapour quenching and 
temperature-dependent collisional cross sections. The 
fourth implementation of a fluorescence model follows 
the calculation from Keilhauer et al. in 2008 |20|. 
Here, 23 wavelengths between 300 and 430 nm are 
considered by applying Eq. (|2]l and (|3]l. The model 
uses a compilation of different measurements |20|, 
II2TI . For the temperature-dependent collisional cross 
sections, the data from AIRFLY |9| are used. These a- 
coefficients are obtained in air, so the same ax is applied 



to NN-collisions and NO-collisions. The temperature- 
dependent collisional cross sections in Eq. (|3]l are written 
as aN:.,.(r) = aS,,,-r"- where a", ^ = aN^^.-293-"- 
is the measured temperature-independent cross section at 
standard experimental conditions of T = 293 K. Cross 
sections for nitrogen - water vapour collisions have been 
measured by two experiments ll22l . ||23| . 

III. Atmosphere-dependent Fluorescence 
Emission 

To study the overall effect of different fluorescence 
models on reconstructed air shower observables, primary 
energy E and position of shower maximum X,„ax, 
it is important to account for only that part of the 
fluorescence spectrum that a detector is sensitive to as 
well as the wavelength dependent attenuation in the 
atmosphere (see for instance Fig. 8d in 1,111). Moreover, 
since the atmospheric parameters P, e and T depend 
on altitude, different fluorescence models will propagate 
differently to E and X^ax if the shower reached its 
maximum high in the atmosphere or close to the ground. 

To include all these effects, we proceeded as follows: 
Proton and iron showers with energies between 10^^ 
and 10^" eV were generated using CONEX |j24l and 
QGSJETII II25I . The fluorescence light was generated 
according to the K08 model including water vapour 
quenching and temperature-dependent collisional cross 
sections. The events were generated with time stamps 
that corresponds to nights with balloon launches, such 
that realistic profiles for P, e and T could be obtained. 
In the following, we will compare the difference in the 
reconstructed E and X^ax values of these simulated 
showers. 

A. Fluorescence Models 

The X,„ax and energy differences for reconstruc- 
tions with different fluorescence models is shown in 
Fig. [T] For this figure, the water vapour quenching and 
temperature-dependent collisional cross sections were 
not switched on in the K08-model, thus this comparison 
is only sensitive to the F(P, T) implementations. As 
explained above, the A07 model is normalised to N04, 
therefore they are not independent and show correspond- 
ingly the smallest differences. 

B. Temperature -dependent collisional cross sections and 
vapour quenching 

The influence of the water vapour quenching (ue) and 
temperature-dependent collisional cross sections {ax) on 
Xmax and E was studied by subsequently switching off 
the effects in the reconstruction using the K08-model. 
As can be seen in the left panel of Fig. |2l ignoring 
(Te- and (TT-effects leads to an underestimation of the 
reconstructed energy by about 5%. Both a^- and ctt- 
dependences affect the shape of the longitudinal profile. 
Since the de -dependence is most important close to 
ground and the ctt -dependence affects mainly higher 
altitudes, the two effects partially compensate (see right 
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Fig. 1. Comparison of 

the influence of different 
ifuorescence models on E 
and Xniax (without vapour 
quenching and temperature- 
dependent collisional cross 
sections). The abbreviations of 
the different fluorescence models 
are defined in Sec. llTl 
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panel of Fig. |2]l leading to only a small X^ax shift of 

< 2 g cm^^. 

Interchanging the water vapour quenching from 1221 
with the independent measurement from |'23'| affects the 
shower observables very little (see solid black dots in 
Fig.©. 

The varying strengths of the a^- and CT-dependences 
at different altitudes can be seen in Fig. [3] Ignoring the 
(It -effect, the energy is misreconstructed up to -7% for 
showers with Xj^ax high up in the atmosphere. Ignoring 
the (Te-dependence, the energy is underestimated also 
up to 7% for showers with X,nax close to ground. 
The position of shower maximum is also affected with 
the largest biases being observed for deep and shallow 
showers. The overall shift of X^ax is strongest for 
showers with a position of shower maximum at about 
3 km a.s.l. with -5 g cm^-^ or for showers with Xmax at 
9 km a.s.l. with 5 g cm^^. It can clearly be seen in the 
right-hand plot of Fig. [3] that the Ue -dependence cancels 
out partly the ctt -dependence concerning X,nax- 

IV. Discussion of Results 

In the fluorescence model K08, all currently known 
effects of the fluorescence light emission are included 
in dependence on varying atmospheric conditions. Run- 
ning this model in combination with actual atmospheric 
profiles, gives a good estimate of the overall mis- 
reconstruction and uncertainties of a standard recon- 
struction. However, it must be stressed that all of the 
models used in this study have a reported uncertainty 
of the absolute fluorescence yield well above 10%. In 
particular, the AIRFLY and AirLight experiments will 
perform an absolute fluorescence yield calibration with 
higher accuracy and results can be expected within one 
year. 

In Fig. m the difference of the reconstruction of 
E and Xmax using the K08 model with all effects 
in combination with actual atmospheric profiles and 
a standard reconstruction with the A07 fluorescence 
model with monthly models can be seen. More or less 
independent of energy, the reconstructed primary energy 
E is higher by about 5% using K08 compared with 
A07 model. The position of shower maximum Xmax is 
nearly unaffected. These results are very similar to the 
comparison of the full K08 model and that without ae- 
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and (iT-dependences. Thus, no additional systematics are 
introduced while changing the fluorescence model apart 
from those obtained by the cTg- and CT-dependences. 

Studying the variation in E and X^ax in dependence 
on the height of the shower maximum, two extreme 
cases can be found: The average shift in E can be up 
to -7% for E and -5 g cm~^ for Xmax for deeply- 
penetrating showers and up to -7% for E and -1-5 g cm^^ 
for Xmax for showers that develop high in the atmo- 
sphere. 

Furthermore, we studied the influence of different 
types of primary particle in terms of proton- and iron- 
induced showers. Comparing the widths of the distribu- 
tion, no difference could be found between proton- and 
iron-induced air showers. 

The change in the atmosphere description from 
monthly models to actual sounding profiles do hardly af- 
fect the reconstructed energy nor the position of shower 
maximum. For E, the difference is well below 1% and 
for Xniax below 2 g cm^^. 

Obviously, the fluctuation of the atmosphere around 
the monthly average atmosphere values adds an ad- 
ditional contribution to the statistical uncertainty of 
the reconstructed energy and X^ax of one shower 
The 'end-to-end' comparison of the A07 model with 
monthly averages to the K08 model with sounding data 
yields RMS{AE/E) e [1.5,3.0]% and RMS(X,nax) e 
[7.2, 8.4] g cm-2(cf. Fig. gjl. 

Finally, the systematic difference in the collisional 
cross section data from two independent measure- 
ments ll22]| . Il23l are negligible. The reconstructed energy 
varies less than 1 % and the position of shower maximum 
about 1 g cm^'^ while interchanging the cross sections. 
Varying the a-coefficients for the temperature-dependent 
collisional cross sections within their given uncertainties, 
yields in less than 1% change in reconstructed energy 
as well. 
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Fig. 2. Comparison of the effect 
of switching off ere and the col- 
lisional cross sections ctt on E 
and Xmax as well as the influ- 
ence of d iffer ent vapour quench- 
ing, a-^ (23) and a^' (2^. 
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Fig. 3. Difference in recon- 
structed energy and Xmax in 
dependence the vertical height 
of the shower maximum (E = 
10^"^ eV). 
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Fig. 4. DiiTerence of the re- 
construction using the full K08 
model to A07. Error bars denote 
the RMS spread. Note that we 
corrected for the 'trivial' yield 
difference. A, from Fig. [T] 
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